Realization of SOC behavior in a dc glow discharge plasma 
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Experimental observations consistent with Self Organized Criticality (SOC) have been obtained 
in the electrostatic floating potential fluctuations of a dc glow discharge plasma. Power spectrum 
exhibits a power law which is compatible with the requirement for SOC systems. Also the estimated 
value of the Hurst exponent (self similarity parameter), H being greater than 0.5, along with an 
algebraic decay of the autocorrelation function, indicate the presence of temporal long-range corre- 
lations, as may be expected from SOC dynamics. This type of observations in our opinion has been 
reported for the first time in a glow discharge system. 



I. INTRODUCTION 

The study of nonequilibrium phenomena of naturally 
occurring [l| , and laboratory plasmas [l], 0| has been an 
active area of research in plasma physics. Within this 
frame work, the Self Organized Criticality (SOC) con- 
cept has been quite rigorously deployed to explain some 
of the turbulent transport observations in magnetically 
confined fusion devices like Tokamaks 0, 0] . The physics 
of glow discharge plasmas [1] in the last two decades, have 
generated a renewed interest due to their importance in 
low temperature plasma applications Being a non- 

linear medium, they have been a good test bed to investi- 
gate various nonlinear phenomena like chaos etc 0, S Q ■ 
In this paper we have attempted to apply, quite success- 
fully the SOC concept for the first time to some of the 
turbulent fluctuations in a glow discharge plasma. 

The outline of the paper is as follows: In Section |TT] 
we describe the experimental setup and observations. In 
Section lnll we present the analysis techniques and results 
of the spectral and statistical methods that have been 
carried to substantiate the SOC behavior. Finally we 
summarize our results in Section ITVl 



II. EXPERIMENTAL SETUP AND RESULTS 

The experiment was carried out in a coaxial cylindrical 
glow dc discharge plasma system with argon as shown in 
Figure [TJ The hollow Stainless Steel (SS) outer cyhn- 
der of 45 mm diameter is the cathode and the SS rod of 
1 mm diameter inside the cathode is the anode, which 
is grounded. The whole system was placed in a vac- 
uum chamber, and evacuated to a base pressure of 10"'^ 
Torr by means of a rotary pump. Argon gas was in- 
troduced using precision needle valve into the chamber. 
Plasma discharges were obtained over a wide range of 
filling pressure and discharge voltage. A Langmuir probe 
made of tungsten of diameter « 0.5 mm and length w 2 
mm was used to measure the electrostatic floating po- 
tential fluctuations in the plasma at about 12.5 mm from 



the anode (Fig. [T]). The fluctuating data of 2500 points 
was recorded at a sampling rate « 1G~^ sec using a Tex- 
tronix oscilloscope, and then transferred to the computer 
for further analysis. 
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FIG. 1: Schematic diagram of the experimental setup of the 
cylindrical dc discharge plasma system with Langmuir probe. 
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FIG. 2: The variation of plasma discharge current (solid line), 
and plasma floating potential (dotted line) with pressure. 

Keeping the discharge voltage at a constant value of 
800 V, a discharge was struck with a very faint glow at 
~ 9 X 10"'^ Torr, and, then as the filling pressure was 
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FIG. 3: Plot of critical pressures with discharge voltages. Cor- 
responding floating potential, and discharge are within brack- 
ets. The first value in the brackets is the floating potential, 
and the second one is the discharge current in (V, mA) unit. 



w —5 V. So the critical pressure divides the scanned re- 
gion into two regions I and II, shown by a vertical line 
(— ) in Figured The plasma density, temperature, and 
electron-electron collision mean free path in region II are 
10"^ - lO^CTO"^, 2-4 eV, and 2.6 x 10^ cm respectively. 
However, in region I, it was almost impossible to obtain 
the I-V characteristics because of their extremely low val- 
ues. Qualitatively looking at discharge current, and glow 
intensity, region II is probably a normal glow discharge 
region, while region I might be the dark or subnormal 
glow discharge region. The critical pressure where the 
transition takes place is not a fixed point, but decreases 
with the discharge voltage as seen in Figure [S] It is also 
seen that the floating potential, and the discharge current 
also decrease with discharge voltage. Figures IH^a) , (b) , 
and (c) are the typical electrostatic fluctuations at pres- 
sures 0.9 X 10~^, 1.2 X 10"^, and 1.5 x 10^^ Torr respec- 
tively in region I, while (d), and (e) are the fluctuations 
at pressures 1.7 x 10^^ and 2.2 x 10^^ Torr respectively 
in region II. 
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FIG. 4: Figure shows the electrostatic fluctuations at pres- 
sures 0.9 x 10~^ (a), 1.2 X 10"^ (b), 1.5 x 10"^ (c), 1.7 x 10"^ 
(d), and 2.2 x 10"^ Torr (e) respectively. 



gradually increased by means of the needle valve, the 
intensity of the glow suddenly became bright at about 
1.6 X 10^^ Torr. The sudden change in the plasma dy- 
namics at the critical pressure led to a change in the float- 
ing potential, and the discharge current. The variation of 
the plasma floating potential, and the plasma discharge 
current with pressure have been shown in Figure [2] It 
shows that, the current (solid line) jumped to a larger 
value by a factor of 15 at 1.6 x 10^^ Torr, and then in- 
creased gradually with pressure. On the other hand. The 
floating potential rapidly fell to a negative value (« —9 
V) up to the same critical pressure (1.6 x 10~^ Torr), 
and then again increased (dotted line in Fig. [J) with fur- 
ther increase in pressure, until it flnally settled down to 



III. ANALYSIS OF SOC BEHAVIOR 

The experimental evidences considered as main ingre- 
dients of SOC are 1//^ {(3 >0) power law (where / 
is the frequency of the fluctuations obtained from Fast 
Fourier Transform) [l^, [HI, [H, , long-range correla- 
tion [Ti] , and nongaussian probability distribution func- 
tion (PDF) 15]. From the power spectral analysis we 
have estimated the f3 from In(Power) versus ln(/). For 
long-range time correlation we estimated the Hurst expo- 
nent H, and the exponent (a) of Autocorrelation function 
(ACF) decay, as described bellow. 

Hurst Exponent-The Rescaled-Range statistics 
{R/S) method was proposed by Hurst and well estab- 
lished by Mandelbrot, and Wallis [lB|. For the time series 
deflned above, the R/S is defined as [H, [l3| the ratio of 
the maximal range of the integrated signal normalized to 
the standard deviation: 



R{n) _ max{0,Wi,W2,...,Wn) -min(0,Wi,W2, 

_ (1) 

Here Wk — xi+X2+X3 + ...+Xk — kX{n), where X, S'^{n), 
and n are respectively the mean, variance, and time lag 
of the signal. The expected value of R/S scales like cn^ 
as n oo, where H is called the Hurst exponent. For 
random data H=0.5, while H > 0.5 for the data with 
long range correlations. H < 0.5 indicates the presence 
of long-range anti-persistency in the data. 

The ACF has been derived as follows: 

Auto-correlation-FoT a time series of length n, X = 
iXj . i = l,2,...n], the ACF function can be written as 



C{t) 



j+T 



X){X,-X) 



(2) 



3 



where X, and r are the mean, and time lag of the time 
series respectively. If there is long-range time dependence 
in the system, then the algebraic decay of the ACF can 
be written as [l^ 

C{t) ^ r-" (3) 

for large t, where < a < 1. 

In order to verify nongaussianity we obtained the PDF 
of the fluctuating data. 
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FIG. 5: In(power) vs. In f plot. The solid line shows the best 
fit. 
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FIG. 6: ln(R/S) vs. In(lag) plot for the electrostatic fluctu- 
ation at the pressure 1.4 x 10~^ Torr. The solid line shows 
best fit. 

Figure [5] shows In (Power) vs. In / of the fluctuations 
in region I from which we estimated the exponent to be 
« 1.6. This agrees well with the numerical (lol . [20j . 
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FIG. 7: R/S as a function of time lag for the electrostatic 
fluctuation of the coherent oscillations at the pressure 1.7 x 
10~^ Torr. For one time period lag, H=l, and for lag more 
than one period, H=0. 
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FIG. 8: ln(R/S) versus In(lag) plot at the pressure 2.2 x 10"^ 
Torr, more than one slope indicates instabilities with many 
modes. 

as well as experimental observations [ll|, [13, in the 
presence of SOC behavior. 

Figure [H] shows a typical plot of ln(R/S) vs. ln(time 
lag) of the fluctuations in region I, for a pressure of 
1.4 X 10"^ The Hurst exponent H is about 0.96 ± .01 
(this indicates long-range time correlations) [2l| . and it 
is almost constant over the entire pressure range in re- 
gion I. On the other hand, for sinusoidal like oscillation 
in region II, the Hurst exponent is 1 for the lag length of 
one period of oscillations, and zero with more than one 
period lag [11] as shown in Figure [71 Also in the same 
region (II) multi slope ln(R/S) vs. In(lag) plot as seen 
in Figure [5] probably due to plasma instabilities of many 
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FIG. 9: In-ln plot of ACF vs. time lag at pressure 1.4 x 10~^ 
Torr. Up to 6 decorrelation times, it shows power law (dotted 
line), and after that it follows exponential decay. 
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FIG. 10: The probability distribution function (PDF) of the 



fluctuation at pressure 1.2 x 10 
the corresponding gaussian fit. 



Torr. and the dotted line is 



same figure. We suspect that there rnight be a slightly 
bimodal distribution similar to Ref . ■ Figure [TUT b) 
shows the Gaussian nature of the fluctuations in region 
II. 

Our results of Hurst exponent, H > 0.5, ACF expo- 
nent, a ~ 0.30, nongaussian PDF, and power spectral 
index (3 ~ 1.60 in the pressure range 9 x lO""^ — 1.6 x 10~^ 
Torr, are consistent with the systems exhibiting SOC like 
behavior. 

Comparison of a, H by ACF, H by R/S, and /? have 
been shown in Table [H for pressures 0.9 x 10~^, 1.2 x 
10^2 1.5 X IQ-^ Torr. 



IV. CONCLUSION 

We have obtained SOC behavior over a finite range 
of neutral pressure of 9 x 10^"^ — 1.6 x 10^^ Torr for a 
fixed discharge voltage. Glow discharges are simple sys- 
tems, but their physics can be quite complicated due to 
the presence of several phenomena like avalanche break- 
down, ionization waves, low frequency ion-acoustic insta- 
bility, double layer, chaos etc. Most of them are highly 
nonlinear processes and hence one requires different tech- 
niques both statistical and spectral to investigate and 
understand their behavior. From our present analysis 
we observe that the plasma dynamics in the region I is 
compatible with self organized criticality, while region II 
is not. This could also imply that plasma transport in 
region I is quite different from region II. Detailed investi- 
gations of the spatio-temporal chaos, and the multifractal 
nature of these fluctuations are in progress and will be 
reported else where. 
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frequencies [1^. The ACF exponent a has been calcu- 
lated for the fluctuations in region I from the In(ACF) vs. 
ln(time lag) plot as shown in Figure^ The ACF (Fig. [5]) 
shows power law up to about 6 times the decorrelation 
time, and after that it follows exponential decay. Average 
value of a is about 0.30. Using the relation H — (2 — a)/2 
[T9| . H calculated from ACF is ~ 0.85, which is close to 
the value of H, calculated using R/S technique. In region 
II no power law decay of ACF has been observed. The 
PDF of the floating potential fluctuations in region I seen 
in Fig. [TUT a) clearly shows a nongaussian nature. Corre- 
sponding best gaussian fit is given by dotted curve in the 



TABLE I: In the following table the decay exponent a of 
the ACF, H from ACF, H using R/S, and the power spectral 
index /3 have been shown in the second, third, fourth, and fifth 
column respectively, for the pressures shown in the column 
one. 

pressure a Hurst Hurst /3 
10-2 (Torr) (ACF) (ACF) (R/S) (PS) 



0.9 


0.38 


0.81 


0.97 


1.85 


1.2 


0.34 


0.83 


0.96 


1.65 


1.5 


0.23 


0.88 


0.95 


1.60 
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